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Abstract 

A  model  of  a  direct  methanol  fuel  cell  anode  is  presented  which  considers  electrocatalysts  to  have  a  distribution  of  over  potential  and 
current  density  in  the  structure.  The  model  is  applicable  to  an  anode  based  on  a  metal  mesh  supported  electrocatalysts  structure.  Methanol 
oxidation  is  described  by  dual  site  mechanisms  involving  adsorbed  CO  and  OH  intermediates.  The  model  is  used  to  predict  the  electrode 
potential-current  density  behaviour  of  the  anode.  The  concentration  of  methanol  is  shown  to  influence  overall  electrode  polarisation 
characteristics  and  critically  the  selection  of  the  mechanism  for  methanol  oxidation  has  a  major  impact  in  this  respect.  The  model  gives 
good  correspondence  with  experimentally  observed  cell  polarisation  behaviour. 
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1.  Introduction 

Proton-exchange  membrane  (PEM)  fuel  cells  are  alterna¬ 
tive  power  sources  for  stationary  and  mobile  applications 
and  electric  vehicles.  Methanol  is  a  liquid  fuel  that  has  sub¬ 
stantial  electroactivity  and  can  be  oxidised  directly  to  carbon 
dioxide  and  water  on  catalytically  active  anodes  in  a  direct 
methanol  fuel  cell  (DMFC); 

CH3OH  +  H20  -+  C02  +  6H+  +  6e“  (1) 

The  direct  methanol  fuel  cell,  based  on  a  solid  polymer  elec¬ 
trolyte  (SPE)  in  the  form  of  a  proton  conducting  membrane, 
has  the  attraction  of  no  liquid  acidic  or  alkaline  electrolyte. 
The  structure  of  the  DMFC  is  a  composite  of  two  porous 
electrocatalytic  electrodes  on  either  side  of  a  solid  polymer 
electrolyte  membrane.  In  the  DMFC,  platinum  alone  is  not 
a  sufficiently  active  methanol  oxidation  electrocatalyst  and 
the  promotion  of  methanol  oxidation  has  been  actively  stud¬ 
ied.  Significant  results  have  been  achieved  with  the  use  of 
binary  catalysts,  notably  Pt-Ru.  With  these  catalysts  the  sec¬ 
ond  metal  forms  a  surface  oxide  used  in  the  potential  range 
for  methanol  oxidation  [1].  The  mechanism  of  methanol  ox¬ 
idation  can  be  represented  as  [2]; 
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CH3OH  -+  COads  +  4H+  +  4e“, 

H20^  OHads  +  H+  +  e-, 

COads  +  20Hads  -+  C02  +  H?0  (2) 

It  is  generally  thought  that  the  rate  controlling  step  is  sur¬ 
face  reaction  between  COads  and  OHads.  It  has  been  pro¬ 
posed  that  methanol  and  hydroxyl  groups  are  adsorbed  on 
different  parts  of  the  surface  on  carbon  supported  platinum. 
However  despite  significant  research  on  methanol  oxidation, 
the  mechanism  is  not  fully  known  and  the  adsorption  of  the 
various  reactive  intermediates  may  involve  a  combination  of 
single  site  and  dual  site  processes  [3], 

Recent  developments  in  electrode  fabrication  techniques 
and  better  cell  designs  have  brought  dramatic  improvements 
in  cell  performance  in  small-scale  DMFCs  [4—14].  An  essen¬ 
tial  condition  for  the  high  performance  of  a  DMFC  is  the  use 
of  relatively  low  methanol  concentrations.  At  concentrations 
higher  than  approximately  2  mol  dm-3,  the  cell  voltage  de¬ 
clines  significantly  due  to  permeation  of  methanol  through 
the  SPE  (Nation®)  membrane,  i.e.  methanol  crossover.  This 
permeation  results  in  a  mixed  potential  at  the  cathode  with 
a  significant  loss  in  oxygen  reduction  performance  and  also 
poor  fuel  utilisation.  Thus  an  important  area  to  improve  the 
DMFC  performance  is  in  polymer  membrane  electrolytes  to 
reduce  methanol  crossover. 

Previous  models  of  the  DMFC  have  been  few.  Verbrugge 
described  a  simple  diffusion  model  of  methanol  through  a 
PEM,  assuming  dilute  solution  theory  [15],  Validation  of 
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the  model  with  experimental  data  showed  that  the  diffusion 
rate  of  methanol  through  the  membrane  was  nearly  as  fast 
as  through  water.  A  second  model,  across  a  PEM  has  been 
used  to  explain  observed  experimental  data  for  a  vapour  feed 
DMFC  [16].  This  model  has  been  extended  to  include  a  one 
dimensional  model  of  the  potential  distribution  and  concen¬ 
tration  distribution  of  methanol  in  the  anode  electrocatalyst 
layer  for  a  vapour  feed  system  [5,17].  The  model  gives  good 
agreement  with  experimental  data  except  under  conditions 
where  mass  transport  becomes  rate  limiting.  Baxter  et  al. 
[18]  have  presented  a  model  of  the  DMFC  anode  which  is 
considered  to  be  a  porous  electrode  consisting  of  an  elec¬ 
tronically  conducting  catalyst  structure  that  is  thinly  coated 
with  an  ion  selective  polymer  electrolyte.  The  pores  of  the 
electrode  are  filled  with  aqueous  methanol  solution  in  which 
all  species  undergo  mass  to  transport.  Mass  transfer  in  the 
anode  is  defined  in  terms  of  a  pseudo-mass-transport  coef¬ 
ficient.  The  model  is  however  not  validated  against  experi¬ 
mental  data  and  does  not  consider  mass  transport  of  species 
in  other  regions  of  the  electrode  assembly. 

A  two-dimensional  mathematical  modelling  of  an 
SPE-DMFC,  both  for  vapour  and  liquid-feed  configura¬ 
tions,  has  been  presented  by  Kulikovsky  et  al.  [19,20] 
employing  Stefan-Maxwell  molecular  diffusion  and  Kniid- 
sen  diffusion  mechanism,  while  neglecting  the  methanol 
crossover  component.  Sundmacher  et  al.  [21]  studied  both 
the  static  and  the  dynamic  response  of  an  SPE-DMFC 
and  showed  that  methanol  crossover  in  the  cell  can  be  re¬ 
duced  by  pulsed  methanol  feed.  Wang  et  al.  [22,23]  have 
presented  a  two-dimensional  and  biphasic,  multicompo¬ 
nent  mathematical  model  for  an  SPE-DMFC.  The  model 
includes  the  crossover  effects  due  to  diffusion,  convection 
and  electro-osmosis  and  the  resultant  mixed  potential  due 
to  methanol  oxidation  at  the  cathode.  Although  the  trans¬ 
port  in  the  porous  media  is  quite  accurately  treated  in  the 
model,  the  catalyst  layer  has  been  considered  to  be  infinites¬ 
imally  thin  excluding  the  polarisation  effects  arising  from 
mass  transport  of  reactants  within  the  layer.  Quite  recently, 
Meyer  and  Newman  [24,25]  have  presented  modelling  and 
data  analysis  of  transport  phenomena  in  an  SPE-DMFC. 
In  contrast  to  most  of  the  earlier  models,  which  employ 
a  simple  Butler-Volmer  (BV)  relationship  for  describing 
the  electrode-kinetics  for  methanol  oxidation  at  the  anode, 
the  model  due  to  Meyer  and  Newman  follows  the  reaction 
mechanism  proposed  by  Gaisteiger  et  al.  [26]. 

In  this  paper  we  present  a  simple  model  of  a  methanol 
oxidation  anode  in  which  a  distribution  in  overpotential  oc¬ 
curs  in  the  structure.  The  model  is  applicable  to  a  new  type 
of  anode  for  the  DMFC  currently  being  investigated  in  our 
laboratories. 


2.  Kinetics  of  methanol  oxidation 

The  following  reaction  mechanism  is  based  on  a  simplifi¬ 
cation  of  a  more  general  model  of  methanol  oxidation  [27]. 


The  mechanism  is  a  reasonable  approximation  in  view  of 
the  great  difficulties  in  determining  kinetic  parameters  for 
the  multi-reaction  paths  of  methanol  oxidation: 


CH3OH4cH3OHads  (3) 

k{ 

CH3OHads4cOad  +  4H+  +  4e"  (4) 

k\ 

H2o4oHads  +  H+  +  e"  (5) 

The  rate  determining  step  (i.e.  slow  process)  is  taken  as 
COads  +  20Hads  h.  C02  +  H20  (6) 


This  will  apply  until  conditions  are  such  that  mass-transport 
limitations  of  methanol  occur,  typically  at  high  current  den¬ 
sities  which  are  not  of  practical  interest. 

We  assume  that  steps  (3)— (5)  are  in  equilibrium.  Step  (5) 
principally  includes  adsorption  on  Ru  sites,  whereas,  steps 
(3)  and  (4)  involve  adsorption  primarily  onto  platinum.  The 
latter  condition  is  valid  as  it  is  known  that  Ru  is  a  poor 
electrocatalyst  for  methanol  oxidation. 

From  equilibrium  of  reactions  (3)-(5)  we  obtain 

k\Cyi  (1  -  0m  -  #C0)  =  &t#M  (7) 


where  9m  and  9qo  refers  to  the  fractional  coverage  of 
methanol  and  CO,  respectively. 

&2#m  =  k\9co  (8) 


where 


k2  =  k®  exp 


(1  -fo)F 

- j 

R 


k\  =  kS]  exp 


E  ' 
-aiF  — 
RT 


noting  that  we  do  not  assume  symmetry  of  electrochemical 
reaction  and 


k3(l  —  0oh)  =  £30oh 
where 


k3  =  k3  o  exp 
k\  —  k\G  exp 


(1  -  03)F — 
H  RT 


E 

-foF— 

RT 


The  activity  of  water  is  assumed  to  be  1 .0. 
Eqs.  (7) — (9)  combine  to  give: 

K2R\CM 


Oco  = 

$OH  = 


(1  +  7Q,  CM(1  +  K2)) 

k3 


1  +  k3 

Now  Kj  ~S>  1  therefore  Eq.  (10)  becomes 
KiK\Cm 

urn  =  - 

1  +  KlK2CM 

where  K\  =  k\/k\,  K3  =  k2/^2,  K3  =  k3/ky 


(9) 


GO) 

(11) 

(12) 
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The  rate  of  reaction  (6)  is  therefore 


J  2 

''6  =  —  =  ^4^6h^CO  =  h 

Of 


k3 


i  +  k3 


K2K1CM 
l  +  K2KiCM 

(13) 


We  define  the  potential  E  in  terms  of  overpotential  rj  and  the 
open-circuit  potential  U0,  i.e. 

E  =  r,+  U0  (14) 

The  reaction  current  density  is  therefore  written  as 
.  _  fci0CMexp(((l  —  j02  +  a2)/RT)Frj ) 

3  1  +  CM&20  exp(((l  -  fi2  +  a2)/RT)Fri) 

where 


k\ 

km  =  6Fk4  — r  exp ( 1  -  +  a2) 

k  "  Arj 


and  k2o  = 


k  10 
k^bF 


FUq 

RT 


and  noting  Kj  1  for  most  practical  fuel  cell  conditions. 

With  the  substitution  ft  =  (1  —  P2  +  <X2)(F/RT)  we  write 
Eq.  (15)  in  a  more  convenient  form  as 

.  6FkmCM^ 

3  1  +  k20CM  j 

or 

.  _ _ Cm _ 

3  [(l/6Ffci0exp03?7))  +  (£2oCm/6FA:io)] 


Fig.  1  shows  experimental  data  for  the  oxidation  of  methanol 
at  a  fuel  cell  anode  at  temperatures  of  60  and  90  °C  [28].  The 
data  clearly  show  that  methanol  oxidation  exhibits  a  limiting 
current  density  which  may  be  due  to  mass  transport  or  to  the 
mechanism  of  methanol  oxidation  in  which  adsorption  is  a 


major  factor.  This  behaviour  is  further  discussed  later.  Plot¬ 
ted  on  the  figure  is  a  Tafel  approximation  of  the  data  based 
on  fitting  the  following  equation  at  low  current  densities: 

j  =  6F£i0Cm  exp  (/}??)  (17) 

Clearly  over  the  lower  range  of  potentials  the  fit  to  the  Tafel 
type  equation  is  reasonable. 

Fig.  2  shows  experimental  data  for  methanol  oxidation 
obtained  at  60  and  90  °C,  and  the  fit  of  the  data  to  Eq.  (16). 
The  fit  has  been  achieved  by  using  the  value  of  the  limiting 
current  density,  as  defined  from  Eq.  (16)  as  j  =  km/ k2o,  to 
determine  k2o  and  then  using  the  data  obtained  for  the  Tafel 
correlation  for  k2o  and  ji.  The  data  at  the  two  temperatures 
has  been  correlated  using  two  values  of  /I  and  two  values 
of  the  constant  k2o .  The  model  shows  reasonable  agreement 
with  experimental  data. 

Fig.  3  shows  the  effect  of  changing  the  value  of  on  the 
correlation  of  the  experimental  data  aimed  at  improving  the 
fit  in  the  range  of  current  densities  >1000 Am-2.  This  is 
done  by  reducing  the  values  of  the  constants  km  and  k2a, 
when  the  value  of  /I  was  reduced.  Although  a  reasonable 
agreement  with  data  is  achieved  at  higher  current  densities, 
in  this  case  the  fit  to  the  data  at  low  current  densities  is  not 
good. 

Fig.  4  shows  the  effect  of  varying  the  ratio  of  k2o/km  °n 
the  anode  polarisation  and  compares  the  resulted  fits  with 
experimental  data  and  the  results  based  on  the  Tafel  approx¬ 
imation  of  Eq.  (16).  The  term  k2o/km  actually  relates  to  the 
rate  coefficient  for  reaction  step  (6)  in  the  mechanism  and 
the  smaller  its  value  the  greater  is  the  rate  of  this  step.  Hence, 
as  expected,  with  a  relatively  fast  reaction  rate,  between  the 
surface  species,  methanol  oxidation  overpotential  behaviour 
becomes  more  favourable,  i.e.  polarisation  is  lower  at  a  given 
current  density,  and  the  limiting  currents  are  higher. 

Fig.  5  shows  the  effect  of  methanol  concentration  on 
anode  polarisation.  At  a  fixed  potential  the  overpotential 


Fig.  1.  Experimental  data  for  the  oxidation  of  methanol  at  a  fuel  cell  anode  at  a  temperature  of:  (A)  90 °C  and  (•)  60 °C  [28]. 
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Fig.  2.  Experimental  data  for  methanol  oxidation  obtained  at  60  and  90°C  [28]  and  the  fit  of  the  data  to  Eq.  (16):  (90°C,  A)  k\ o  =  8.64  x  10  l0, 
k2 0  =  5.83  x  lO”8,  P  =  24.91;  (60°C,  O)  *10  =  8.64  x  10"10,  k20  =  1.21  x  1(T7,  p  =  21.84.  CM  =  lOOOmolm-3. 
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Fig.  3.  The  effect  of  p  on  the  correlation  of  experimental  data  using  model  equation  (16)  for  methanol  oxidation:  (90°C,  A)  kio  =  1.22  x  10  8, 
P=  17.92,  k20  =  7.68  x  10"7;  (60°C,  •)  A10  =  1.22  x  1(T8,  p  =  14.95,  k20  =  1.45  x  10-6. 


decreases  as  the  methanol  concentration  is  increased.  The 
data  exhibits  the  characteristic  that,  at  high  potentials,  a 
common  limiting  current  is  achieved  at  all  concentrations 
of  methanol.  This  behaviour  in  contrary  to  that  reported  in 
a  number  of  studies  of  methanol  oxidation  and  may  be  a 
result  of  the  influence  of  mass-transport  limitations  in  prac¬ 
tical  fuel  cell  anodes  or  to  limitations  in  the  model  above. 
We  address  this  latter  point  first. 

2.1.  The  Gasteiger  mechanism 

An  alternative  mechanism  for  methanol  oxidation  [26]  is 
similar  to  that  proposed  above 


CH3OH  %  CH3OHads  (18) 

*J 

CH3OHads4cOad  +  4H+  +  4e"  (19) 

k\ 

H2C)40Hads  +  H+  +  e"  (20) 

k\ 

The  rate  determining  step  (i.e.  slow  process)  is  taken  as 
COads  +  20Hads  C02  +  H20  (21) 

except,  when  the  concentration  of  methanol  is  small. 
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Fig.  4.  The  effect  of  varying  the  ratio  of  ^20/^10-  (A)  100,  (A)  75,  (O) 
50,  (#)  25,  (+)  10.  Dotted  line:  experimental  data;  continuous  line:  Tafel 
equation. 
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Fig.  7.  The  effect  of  methanol  concentration  on  the  methanol  oxidation 
polarisation  characteristics:  (O)  0.125,  (A)  0.25,  (O)  0.5,  (#)  1.0,  (♦) 
1.5,  (A)  2.0.  90  °C.  kw  =  1.22E-05,  p  =  17.91819,  k2 0  =  0.768343. 
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Fig.  5.  The  effect  of  methanol  concentration  on  the  polarisation  curve 
for  the  methanol  fuel  cell  anode:  (O)  125,  (#)  250,  (A)  500,  (O) 
1500,  (A)  2000.  *2o/*io  =  70.  90°C.  kl0  =  1.22  x  10~8,  p  =  17.92, 
k20  =  7.68  x  10"7. 


The  model  assumes  that  reaction  steps  (19)  and  (20)  are 
in  equilibrium  and  produces  the  overall  expression: 

.  _  6Fk\oCM  exp(/fy) 

3  CM  +  k2  0  exp(0?7) 

where  the  terms  k\o  and  kij)  are  different  group  terms  to 
those  in  Eq.  (16). 

Fig.  6  shows  the  fit  of  model  equation  (22)  to  experimental 
data  at  temperatures  of  60  and  90  °C.  With  the  data  available, 
this  model  also  gives  a  good  prediction  of  methanol  oxida¬ 
tion  polarisation  characteristics.  However  Eq.  (22)  predicts 
(Fig.  7)  different  behaviour  to  Eq.  (15),  i.e.  at  low  overpo¬ 
tential,  the  current  density  at  a  fixed  overpotential  is  inde¬ 
pendent  of  Cm  and  at  high  overpotential  j  is  proportional  to 
Cm-  There  is  a  body  of  literature  data  that  agrees  partly  with 
the  latter  characteristics.  However  there  is  less  evidence,  in 
terms  of  reported  data,  that  is  in  accord  with  the  invariance 
of  methanol  oxidation  currents  with  methanol  concentration, 
at  low  overpotentials. 
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Fig.  6.  The  comparison  of  experimental  methanol  oxidation  data  with  the  model  based  on  the  mechanism  of  Gasteiger:  (90  °C,  A)  k io  =  1.22  x  10  \ 
P  =  17.92,  k20  =  0.77;  (60 °C,  •)  *i0  =  1.22  x  10~5,  p  =  14.95,  k20  =  1.45. 
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2.2.  Implications  of  diffusion  mass  transport 


The  basic  difference  between  the  two  models  of  the 
mechanism  for  methanol  oxidation  presented  above  is  that 
in  one  case  the  adsorption  of  methanol  species  is  taken 
as  an  equilibrium  process  (Frelink)  and  in  the  other  case 
it  is  not  (Gasteiger).  The  consequence  of  this  difference 
is  that  a  limit  to  the  current  density  for  methanol  oxida¬ 
tion  arises  which  in  one  case  is  independent  of  methanol 
concentration  and  in  the  other  case  it  is  approximately 
proportional  to  the  methanol  concentration.  In  previous  ex¬ 
perimental  work  [7,29]  with  fuel  cell  electrodes,  limiting 
currents  have  been  produced  which  increase  with  an  in¬ 
crease  in  methanol  concentration.  These  limiting  currents 
have  been  interpreted  as  either  due  to  mass-transport  limi¬ 
tations  in  the  fuel  cell  electrode  [29]  or  as  a  consequence  of 
the  mechanism  of  methanol  oxidation  [2].  However,  with 
much  of  the  published  data  for  fuel  cell  electrodes,  based 
on  carbon  supported  catalysts,  it  is  not  readily  possible  to 
differentiate  the  major  factor  which  causes  a  limiting  cur¬ 
rent.  A  factor  in  the  DMFC  is  that  of  the  gas  evolution  in 
the  porous  structure  which  will  have  a  controlling  influence 
on  diffusion  and  convective  mass  transport.  Essentially  as 
the  methanol  concentration  is  increased,  then  any  increase 
in  limiting  current  may  be  a  consequence  of  the  greater 
rate  of  carbon  dioxide  evolution.  Furthermore  it  is  also 
possible  that  both  kinetic  (mechanistic)  and  mass-transport 
limitations  may  occur  at  similar  overpotentials  or  cur¬ 
rent  densities.  For  example,  if  we  consider  the  limit  of 
the  Gasteiger  model  at  higher  overpotentials,  then  we  see 
that 


6F7:ioCm 

&20 


(23) 


where  Cm  is  the  concentration  of  methanol  at  the  catalyst 
surface. 


For  a  simple  mass-transport  model,  based  on  a  mass- 
transport  coefficient  k\  ,  the  current  density  can  be  expressed 
as 

j  =  6FkL(CMb  -  Cm)  (24) 

where  CMb  is  the  bulk  concentration  of  methanol. 
Combining  these  two  equations  we  obtain 

- 6  FkLCMb - 

J  (1/*l)  +  (*20/*10) 

The  limiting  current  is  proportional  to  the  methanol  concen¬ 
tration  and  its  magnitude  depends  upon  the  combined  values 
of  the  kinetic  and  mass-transport  coefficients. 

To  determine  whether  the  limiting  current  is  a  result  of 
the  reaction  mechanism,  requires  a  study  of  methanol  ox¬ 
idation  under  well  defined  and  controllable  mass-transport 
conditions,  for  example  using  a  rotating  disc  electrode. 

Furthermore,  the  majority  of  methanol  fuel  cell  anodes 
are  high  surface  area  materials  which  will  exhibit  a  current 
distribution  normal  to  the  direction  of  current  flow.  This  fac¬ 
tor  will  also  have  an  influence  on  the  methanol  anode  polar¬ 
isation  behaviour  as  is  the  subject  of  the  following  section. 


3.  Current  distribution  model 

We  adopt  a  model  for  the  catalyst  structure  in  which  the 
supported  catalyst  is  electronically  in  contact  with  the  cur¬ 
rent  collector  e.g.  carbon  and  ionically  in  contact  with  the 
membrane  through  a  Nation  bonding  layer.  The  concentra¬ 
tion  of  methanol  within  the  pores  is  assumed  uniform. 

The  region  from  x  =  0  to  x  =  1  is  the  catalyst  region 
in  which  the  concentration  of  reactant  is  assumed  constant 
(Fig.  8).  This  situation  applies  to  a  system  with  negligi¬ 
ble  mass-transport  resistance  in  the  pores  or  when  an  ap¬ 
proximate  well  mixed  condition  applies.  This  applies  to  an- 


Fig.  8.  Schematic  model  of  methanol  fuel  cell  anode. 
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odes  which  are  based  on  Ti  mesh  supported  electrocatalysts 
[30], 

Applying  Ohm’s  law  to  proton  and  electron  motion  in 
the  carbon/catalyst  region  and  ionomer  region  covering  the 
catalyst  gives 

dcpm 

— a - 

dx 


=  J 


(26) 


-K—  = 
dx 


J  s 


(27) 


where  j  and  js  denote  proton  and  electronic  current  den¬ 
sity,  respectively;  <f>m  and  (j)s  the  electrical  potential  in  the 
ionomer  and  carbon  phase;  and  a  and  x  the  corresponding 
effective  conductivities. 

The  proton  and  electron  current  densities  are  related  to 
the  respective  molar  fluxes  according  to 

—  =  (28) 

dx  dx 

Integrating  over  the  catalyst  layer  gives 

j  +  h  =  j\  (29) 

where  j\  is  the  current  density  loading  for  the  cell,  noting 
that  at  x  =  0,  the  membrane  is  electronically  insulating,  i.e. 
j  U=0  =  0  and  thus  j\  =  j\x=0 

In  the  catalyst  region  the  over  potential, driving  reaction 
is  defined  as 

b  =  -(</>s  -  </>m)  (30) 

The  change  in  over  potential,  from  Eqs.  (26)  and  (27),  is 

dn  ( j  /s  \ 

l  j  j s  ,  (31) 


(32) 


dx  \a  K  , 

Substituting  Eq.  (29)  into  (31)  gives 


d  i] 
dx 


1 

~K 


1 


A 

K 


In  many  cases  e.g.  for  carbon  supports,  the  electrical  con¬ 
ductivity  in  the  ionomer  layer  is  much  less  than  that  in  the 
carbon  layer  and  the  terms  1  IK  can  be  ignored  in  Eq.  (32), 
except  at  low  current  densities.  However  other  lower  con¬ 
ductivity  supports  are  being  considered  in  certain  fuel  cell 
(e.g.  conducting  polymers)  and  thus  we  do  not  adopt  this 
assumption. 

3.1.  Model  solution 

In  the  first  instance  we  seek  a  solution  to  the  current 
distribution  problem  assuming  kinetics  for  the  reaction  given 
by  Eq.  (16).  From  this  we  will  establish  the  overall  electrode 
potential  current  density  relationship  for  the  electrode. 

In  the  catalyst  region,  the  change  in  ionic  current  density, 
i.e.  due  to  proton  conductivity  is  given  by 

d  j 


dx 


=  ~aj 


(33) 


where,  a ,  is  a  specific  reaction  surface  area  (m  [). 


Differentiating  Eq.  (16)  we  obtain 
d  j  d?7 


dx  dx 


/>i- 


—tf 

6Fk10 


(34) 


Combining  Eqs.  (33)  and  (34)  we  obtain 


d  ~j  d  j 

XT  =  =  aP 

axz  ax 


=af[  i 


kio  =  6Fkio 
and  hence 

d 2J_  =  _ 

dx2 


•  k20  .2 

J~^J 
k  to  . 

2  _  %)  .3 

ho 


ci/3 

~  -j^jt 


M 

K 


■  k20  .2 

J~ir] 

k\ o  - 


(35) 


1  \  d  j  /  .  k2 o  ,2 
a  J  dx  V  klQ  ^ 


Pje 

K 


j  k20  .  d J 

kio  J  dx 


(36) 


We  now  introduce  dimensionless  variables 


x 

le' 1  =  T 

and 

also 

•* 

e  =  Jo 


—  H -  =7 


Pt-jl 

K 


1 


1  +  (K/o) 

Then  Eq.  (36)  becomes 

j 


d 2j*  _  .»  d j* 
dz2  J  dz 


*2  dj* 


■  e 


r 


!*  dz 

*  r\  i*  _  ti  /r)\  A  ( 


\  _  n  d  r 


r 


f-  (37) 

dz 


noting  that  j*  dj*  —  (1/2)  d (j*)  gives 

d  2r  ,  Id  o'*)2  1  d(/*)3  edj* 


*2 


dz2  2  dz 

Integration  gives 

;*2 


3/*  dz 


dz 


g  dj" 
2 i*  dz 


=  0 


(38) 


d r 

dz 


J 


0  +  f) 


;*3 


J 

3  i* 


ej 


—  _  j* 


j*  3 

Jo 


3  J] 


•sp  I 

g/o  + 


Jo2 


e 

it 


1  +  ^1  =  ^ 


(39) 


with 


f  =  pnlK  +  l-)j 


and 
J*  = 


d f 

dz 


=  aiji 
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where  jg  is  the  current  density  at  the  membrane  catalyst 
interface  which  is  a  function  of  the  electrode  overpotential, 
i.e.  J*  =  al6F(Ckl0e^/(l  +  Ck2 0e^")) 


Re-arranging  Eq.  (39)  gives  the  current  distribution: 


d r 

d  zr 


e  j*1 

1  4“  —  -}-  - 

a£  3  ati* 


=  J  - 


*2 


0  +  £) 


(40) 


Integration  gives 


1 


tan 


,y/-Jd 

a  =  (!  +  Jl)  'd  =  \  0  +  T* )  ’  b 

For  4  =  0  Eq.  (41)  becomes 


Jo 


fit 


:  tan 


J 


~(b/d)j*+  1 


(1  /f^Jd)(d+  (bJ/d))_ 


(42) 


when  e  =  0,  i.e.  a  high  electronic  conductivity  catalyst  layer 
equation  (32)  gives 

1  (  1  -  (2bj*/\) 


Jo  ~ 


:  tan 


yTTJTc  \(1  +  (2fec/l))(l/V— c/2) 


3.2.  Current  distribution  model  for  the  Gasteiger 
mechanism 


By  following  the  same  procedure  as  above  for  developing 
the  current  distribution  model  we  can  generate  the  differen¬ 
tial  equation  defining  the  current  distribution  in  the  anode 
catalyst  as 


d2./  d  j 

TT  =  ~fl  —  =  aP 
dxz  ax 

(  1  1 
-aP\K  +  7 


.  *20  .2 
J-J-J 

*10  J  L 
■2  _  ^20  3 
*10 


1 


a/3 


Je 

K 


*20  .2 
*10  . 


(25') 


where 

*io  =  6ECm*io 

Solution  of  this  equation  using  similar  procedures  as 
above  gives  the  current  distribution  as 

.*  =  1  fan  /  1  ~  (bf0/d) 

Jo  ^dfc  “  V(fl  +  (bc/d))(l/V^cd) 
where 


/,  e  \ 

1  /  e  \ 

,  d  = 

V  at) 

2  V  i*J 

1 

b  — - , 

c  —  J 

3  ati* 

4.  Results  of  current  distribution  model 

The  treatment  of  the  experimental  data  in  Section  2  is 
based  on  an  assumption  that  the  electrode  was  pseudo-two- 
dimensional  and  was  not  porous.  In  practice  the  electrode 
was  made  as  a  Nafion  bound  carbon  supported  structure  sev¬ 
eral  microns  in  thickness  and  exhibits  a  current  distribution 
normal  to  the  flow  of  current.  In  this  section,  we  investigate 
the  simple  current  distribution  model  presented  in  Section  3 
and  its  influence  on  the  overall  anode  polarisation  charac¬ 
teristics. 

Fig.  9  shows  typical  dimensionless  current  density  distri¬ 
bution  in  the  anode  catalyst.  The  behaviour  reflects  the  ex¬ 
pected  decrease  in  local  current  density  from  the  membrane 
surface  into  the  electrocatalyst  structure  associated  the  grad¬ 
ual  decrease  in  ionic  current. 

Fig.  10  shows  typical  dimensionless  polarisation  charac¬ 
teristics  of  the  anode;  i.e.  dimensionless  current  density  as 
a  function  of,  J*,  the  volumetric  current  density. 

Fig.  11  shows  the  effect  of  methanol  concentration  on 
the  overall  anode  polarisation  characteristics  for  the  Frelink 


Fig.  9.  Anode  current  density  distribution:  (•)  J *  =  15,  (A)  J*  =  5,  Cm  =  lOOOmolm  3,  Atjo  =  8.64  x  10  10,  /i  =  24.9,  *>0 
i  =  IE— 5,  a  =  3.4,  K  =  40. 


375, 
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Fig.  10.  Dimensionless  anode  polarisation  behaviour:  (A)  K  =  40,  (•)  K  =  4000,  a  =  130000,  l  =  IE— 5,  a  =  3.4.  k\o  =  1.22E— 05,  f  =  17.91819, 
*20  =  0.768343. 


Fig.  1 1.  The  effect  of  methanol  concentration  on  anode  polarisation  characteristics:  *io  =  8.64E— 10,  fl  =  24.90,  *20  =  5.83E— 08,  a  =  375,  l  =  IE— 5, 
or  =  3.4,  K  =  40.  (♦)  CM  =  2500,  (•)  CM  =  500,  (A)  CM  =  1000,  (O)  CM  =  250. 


model  (Eq.  (16)).  The  values  of  parameters  selected  (thick¬ 
ness,  conductivity)  are  those  that  are  in  general  accord  with 
the  actual  fuel  cell  anode  although  the  exact  values  are  not 
known.  The  model  generally  gives  reasonable  agreement 
with  experimental  data  by  adjusting  the  parameter  *10.  The 
values  of  the  parameter  *10  are  much  lower  than  those  used 
to  model  data,  for  a  pseudo-two-dimensional  electrode,  in 
Section  2  and  is  essentially  due  to  the  much  higher  surface 
area  for  oxidation  in  the  porous  electrocatalyst.  The  data 
show  that  increasing  methanol  concentration  reduces  the  ex¬ 
tent  of  anode  polarisation  and  that  regardless  of  the  methanol 


concentration  one  “limiting”  current  density  is  approached. 
The  behaviour  at  low  polarisation  is  approximately  in  agree¬ 
ment  of  observations  that  the  methanol  oxidation  is  a  half 
order  reaction  in  methanol  concentration,  i.e.  at  a  fixed  over¬ 
potential,  an  order  of  magnitude  increase  in  concentration 
produces  an  approximate  3  to  4  increase  in  current. 

Fig.  12  shows  the  effect  of  methanol  concentration  on  the 
overall  anode  polarisation  characteristics  for  the  Gasteiger 
model  (32).  This  model  with  appropriate  parameters  gives 
polarisation  behaviour  that  is  in  accord  with  experimental 
data.  The  model  exhibits  the  trends  that  at  low  current  den- 
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Fig.  12.  The  effect  of  methanol  concentration  on  anode  polarisation  characteristics,  kio  =  8.46E— 10,  ft  =  24.90,  &20  =  5.83E— 8,  a  =  375,  £  =  IE— 5, 
a  =  3.4,  K  =  40.  (O)  CM  =  2500,  (♦)  CM  =  1500,  (•)  CM  =  500,  (A)  CM  =  1000. 


sities  the  anode  polarisation  does  not  vary  with  methanol 
concentration  and  that  at  higher  anode  polarisation  limit¬ 
ing  current  densities  are  approached  which  depend  upon  the 
concentration  of  methanol.  The  former  observation  is  not 
consistent  with  reported  data  for  methanol  oxidation  by  a 
number  of  researchers  [2,7]. 


5.  Conclusions 

A  model  of  the  current  distribution  in  a  methanol  oxi¬ 
dation  fuel  cell  anode  is  presented  which  can  be  used  in 
conjunction  with  a  kinetic  model  of  the  oxidation  to  pre¬ 
dict  anode  fuel  cell  polarisation  behaviour.  The  ability  of  the 
model  to  accurately  predict  the  polarisation  behaviour  re¬ 
lies  on  the  availability  of  accurate  parameters  for  the  anode 
structure  and  the  kinetic  parameters.  In  addition  the  suitabil¬ 
ity  of  the  model  requires  confirmation  of  the  mechanism  for 
methanol  oxidation,  albeit  an  approximation,  so  that  a  suit¬ 
able  kinetic  equation  can  be  used.  The  latter  equation  should 
be  in  agreement  with  the  effect  of  methanol  concentration 
on  the  kinetics  and  mass  transfer  characteristics  during  ox¬ 
idation.  The  latter  aspect  is  the  subject  of  ongoing  work  in 
the  laboratories  at  Newcastle. 
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